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fSj ■ Abstract 

> ■ 

Cfj , A new evaluation of the hadronic vacuum polarization contribution to the muon magnetic mo- 

ment is presented. We take into account the reanalysis of the low-energy e + e~ annihilation 
cross section into hadrons by the CMD-2 Collaboration. The agreement between e + e~ and r 
spectral functions in the tttt channel is found to be much improved. Nevertheless, significant 
discrepancies remain in the center-of-mass energy range between 0.85 and 1.0 GeV, so that we 
refrain from averaging the two data sets. The values found for the lowest-order hadronic vacuum 
polarization contributions are 



9*' hadLO f (696.3 ± 6.2 cxp ± 3.6 rad ) 10- 10 [e+e" -based] , 



a 



■ a had,LO _ J u """exp 

" " \ (711.0 ± 5.0 cxp ± 0.8 rad ± 2.8 SU(2) ) 10- 10 [r-based] , 

where the errors have been separated according to their sources: experimental, missing radiative 
corrections in e + e~ data, and isospin breaking. The corresponding Standard Model predictions 
for the muon magnetic anomaly read 

| (11 659 180.9 ± 7.2 had ± 3.5 LBL ± 0.4 QED +ew) 10~ 10 [e+e" -based] , 
a "~\ (11 659 195.6 ± 5.8 had ± 3.5 LBL ± 0.4 QED +ew) KT 10 [r-based] , 

where the errors account for the hadronic, light-by-light (LBL) scattering and electroweak con- 
tributions. The deviations from the measurement at BNL are found to be (22.1 ± 7.2 ± 3.5 ± 
8.0) 10~ 10 (1.9 er) and (7.4 ± 5.8 ± 3.5 ± 8.0) 10" 10 (0.7 a) for the e+e"- and r-based estimates, 
respectively, where the second error is from the LBL contribution and the third one from the 
BNL measurement. 
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1 Introduction 



Hadronic vacuum polarization in the photon propagator plays an important role in the 
precision tests of the Standard Model. This is the case for the muon anomalous magnetic 
moment = (g M — 2)/2 where the hadronic vacuum polarization component, computed 
from experimentally determined spectral functions, is the leading contributor to the un- 
certainty of the theoretical prediction. 

Spectral functions are obtained from the cross sections for e + e~ annihilation into 
hadrons. The accuracy of the calculations has therefore followed the progress in the 
quality of the corresponding data [1]. Because the latter was not always suitable, it was 
deemed necessary to resort to other sources of information. One such possibility was the 
use of the vector spectral functions [2] derived from the study of hadronic r decays [3] 
for the energy range less than m T c 2 ~ 1.8 GeV. Also, it was demonstrated that pertur- 
bative QCD could be applied to energy scales as low as 1-2 GeV [4], thus offering a way 
to replace poor e + e" data in some energy regions by a reliable and precise theoretical 
prescription [5-9]. 

A complete analysis including all available experimental data was presented in Ref. [10], 
taking advantage of the new precise results in the 7ttt channel from the CMD-2 exper- 
iment [11] and from the ALEPH analysis of r decays [12], and benefiting from a more 
complete treatment of isospin-breaking corrections [13,14]. In addition to these major 
updates, the contributions of the many exclusive channels up to 2 GeV center-of-mass 
energy were completely revisited. It was found that the e + e~ and the isospin-breaking 
corrected r spectral functions were not consistent within their respective uncertainties, 
thus leading to inconsistent predictions for the lowest-order hadronic contribution to the 
muon magnetic anomaly: 



The quoted uncertainties are experimental, missing radiative corrections to some e + e~ 
data, and isospin breaking. The leading contribution to the discrepancy originated in 
the mi channel with a difference of (—21.2 ± 6.4 cxp ± 2.4 rad ± 2.6su(2) (±7.3 to tai)) 10~ 10 . 
The estimate based on e + e~ data has been confirmed by another analysis using the same 
input data [15]. When compared to the world average of the muon magnetic anomaly 
measurements, 



which is dominated by the 2002 BNL result using positive muons [16], the respective 
e+e~-based and r-based predictions disagreed at the 3.0 and 0.9 a level, respectively, 
when adding experimental and theoretical errors in quadrature. 

The purpose of this letter is to update our analysis [10] in light of the following 
developments. 




[e + e —based] , 
[r— based] , 





(2) 
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• The CMD-2 Collaboration at Novosibirsk discovered that part of the radiative treat- 
ment was incorrectly applied to the data. A complete reanalysis has been carried out 
and presented for publication [17]. As the CMD-2 data dominate the e + e~-based 
prediction (1), the changes produce a significant effect in the final result. Recently 
available results from the SND Collaboration are also included. 

• No significant change occurred for the r-based prediction. The only relevant fact is 
a new result [18] for the branching ratio of the r~ — > v T h~7r mode (h~ stands for 
a charged pion or kaon). 



2 Muon Magnetic Anomaly 

It is convenient to separate the Standard Model (SM) prediction for the anomalous mag- 
netic moment of the muon into its different contributions, 

Q SM = a QED + fl had + ^ (3) 

with 

had had,LO , had, HO , had,LBL / a\ 

and where a® ED = (11658 470.6 ± 0.3) 10~ 10 is the pure electromagnetic contribution 
(see [19,20] and references therein 2 ), a£ ad ' LO is the lowest-order contribution from hadronic 
vacuum polarization, a^ ad > HO = (—10.0 ± 0.6) 10~ 10 is the corresponding higher-order 
part [23,2], and a™ eak = (15.4 ± 0.1 ± 0.2) 10~ 10 , where the first error is the hadronic 
uncertainty and the second is due to the Higgs mass range, accounts for corrections due 
to exchange of the weakly interacting bosons up to two loops [24]. For the light-by- light 
(LBL) scattering part we add the values for the pion-pole contribution [25-27] and the 
other terms [26,27] to obtain a h ^ d ' LBL = (8.6 ± 3.5) 10~ 10 . 

Owing to the analyticity of the vacuum polarization correlator, the contribution of the 
hadronic vacuum polarization to a M can be calculated via the dispersion integral [28] 

ar .™ = ?mj ds m Rish (5) 

where K(s) is the QED kernel [29] , 

K{s) = x 2 (l-y) + ( 1 + ^) 2 ( 1 + ^) + + + [r^ 2 lnx > ( 6 ) 

with x = (1 - /^/(l + and ^ = (1 - Aml/s) 1 ' 2 . In Eq. (5), R(s) = RW(s) 
denotes the ratio of the 'bare' cross section for e + e~ annihilation into hadrons to the 

2 Some adjustment was recently made concerning the fourth-order contribution from the leptonic light- 
by-light scattering, mostly affecting the QED prediction for a e and through it the value of a [21,22]. The 
resulting change in a^ ED is within the quoted uncertainty of 0.3 10~ 10 and has not been included in the 
present analysis. 
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pointlike muon-pair cross section. The 'bare' cross section is denned as the measured 
cross section, corrected for initial-state radiation, electron-vertex loop contributions and 
vacuum polarization effects in the photon propagator (note that photon radiation in the 
final state (FSR) is included in the 'bare' cross section). The reason for using the 'bare' 
(i.e. lowest order) cross section is that a full treatment of higher orders is anyhow needed 
at the level of a^, so that the use of 'dressed' cross sections would entail the risk of double- 
counting some of the higher-order contributions. 

The function K(s) decreases monotonically with increasing s. It gives a strong weight 
to the low energy part of the integral (5). About 91% of the total contribution to a^ ad ' LO 
is accumulated at center-of-mass energies y/s below 1.8 GeV and 73% of a ^ ad ' LO is covered 
by the two-pion final state which is dominated by the p(770) resonance. 

3 Changes to the Input Data 

3.1 e + e~ Annihilation Data 

The CMD-2 data, published in 2002 for the tttt channel [11], have been completely reana- 
lyzed [17] following the discovery of an incorrect implementation of radiative corrections 
in the analysis program. Overall, the pion-pair cross section increased by 2.1% to 3.8% in 
the measured energy range (cf. Fig. 1), well above the previously quoted total systematic 
uncertainty of 0.6%. Specifically, the leptonic vacuum polarization contribution in the 
t-channel had been inadvertently left out in the calculation of the Bhabha cross section. 
This effect produced a bias in the luminosity determination, varying from 2.2% to 2.7% in 
the 0.60-0.95 GeV energy range. The problem consequently affected the measured cross 
sections for all hadronic channels. Another problem was found in the radiative corrections 
for the muon-pair process, ranging from 1.2% to 1.4% in the same region. 
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Figure 1: Relative change in the e + e — > tt + tt cross section of the revised CMD-2 anal- 
ysis [17] with respect to the one previously published [11]. 
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A more refined treatment of hadronic vacuum polarization was performed, with changes 
not exceeding 0.2% for most data points. The effects in the Bhabha- and muon-pair chan- 
nels also affected the event separation. In the CMD-2 analysis, Bhabha events are well 
identified using the electron calorimeter signature while pions are not separated from 
muons. The numbers of electron-, muon- and pion-pair events are based on a likelihood 
method keeping the ratio of muons to electrons fixed through the corresponding QED 
cross sections. Thus, the corrections to these cross sections had an effect on the event 
separation and the measured ratio of pion pairs to electron and muon pairs changed by 
typically 0.7%. 

The correction of the bias in the luminosity determination increases all hadronic cross 
sections published by CMD-2. The changes are 2.4% and 2.7% on the u and resonance 
cross sections, respectively. They are not yet available for the energy range above the 
0. Instead we use an estimated correction of 1.7%, which insignificantly affects the con- 
tribution of the multihadron processes between 1.05 and 1.40 GeV. All these luminosity 
corrections have been applied to the present analysis. Also, the CMD-2 Collaboration 
now provides hadronic vacuum polarization-corrected ui and <f> cross sections so that we 
do not apply this correction anymore. 

Newly published data by SND on the u resonance [30] and the 2ir + 2ir~ as well as 
7T + 7r~27r° modes [31] (unchanged cross sections for the latter two, but reduced systemat- 
ics with respect to previous publications) have been included in this update. 

We refer to our previous analysis [10] for a detailed discussion of radiative corrections, 
in particular the effect of final-state radiation by the charged hadrons. Also given therein 
is a compilation of all input data used to evaluate the integral (5). 

3.2 Data from Hadronic r Decays 

The only update here relates to the normalization of the spectral function in the 7T7t 
channel. New results have been presented by the L3 Collaboration on branching ratios 
for hadronic r decays [18]. Their value for the t' -> u T h'n° mode, (25.89±0.16±0.10)%, 
gives, after correcting for the K~n° contribution [32,33], a result of (25.44±0.16±0.10)% 
for the 7r~7T° mode, in agreement with the previous measurements [12,34,35], yielding the 
world average (25.46 ± 0.10)%, which is used in the present analysis. 

To use the r spectral functions in the vacuum polarization dispersion integral, a value 
for the CKM matrix element \V u d\ is necessary. In the previous analysis, we used the 
average of two determinations [36], |T4<i| = 0.9734 ± 0.0008 from (5 decays and \V U( i\ = 
0.9756 ± 0.0006 from decays and CKM unitarity, which are not consistent. The final 
error was scaled up correspondingly. The determination of V us from hyperon decays [37] 
is in fact more consistent with (3 decays, yielding from unitarity \V ud \ = 0.9744 ± 0.0006. 
New information is expected from recent and neutron decay experiments. For the 
moment we keep our previous average, \V u d\ = 0.9748 ± 0.0010, since the enlarged error 
covers the range of measured values. The V u d uncertainty corresponds to a shift of the 
r-based aj^ ad estimate of 1.1 10~ 10 , which is small compared to the total uncertainty of 
5.8 10- 10 . 
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Figure 2: Relative comparison of the tt + tt spectral functions from e + e and isospin- 
breaking corrected r data, expressed as a ratio to the r spectral function. The band shows 
the uncertainty on the latter. The e + e~ data are from CMD-2 [17], CMD [40], OLYA [40, 
41] and DM1 [42]. The lower figure emphasizes the p peak region. 
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4 Comparison of e + e and r Spectral Functions 

The new e + e~ and the isospin-breaking corrected r spectral functions can be directly com- 
pared for the 7T7T final state. The r spectral function is obtained by averaging ALEPH [3], 
CLEO [38] and OPAL [39] results [10]. The e + e~ data are plotted as a point-by-point 
ratio to the r spectral function in Fig. 2, in a wide energy range (upper plot) and in the 
region around the p peak (lower plot). The central band in Fig.2 gives the quadratic sum 
of the statistical and systematic errors of the r spectral function obtained by combining 
all r data. The e + e~ data have moved closer to the r results: they are now consistent 
below and around the peak, while, albeit reduced, the discrepancy persists for energies 
larger than 0.85 GeV. 

A convenient way to assess the compatibility between e + e~ and r spectral functions 
proceeds with the evaluation of r decay fractions using the relevant e + e~ spectral func- 
tions as input. All the isospin-breaking corrections detailed in Ref. [10] are included. This 
procedure provides a quantitative comparison using a single number. The weighting of 
the spectral function is however different from the vacuum polarization kernels. Using 
the branching fraction £>(r~ — > v T e~v e ) = (17.810 ± 0.039)%, obtained assuming leptonic 
universality in the charged weak current [12], the prediction for the nn channel is 

B C vc(r- -> u T n-7r°) = (24.52 ± 0.26 cxp ± 0.11 rad ± 0.12 su(2) )% , (7) 

where the errors quoted are split into uncertainties from the experimental input (the e + e~ 
annihilation cross sections) and the numerical integration procedure, the missing radiative 
corrections applied to the relevant e + e~ data, and the isospin-breaking corrections when 
relating r and e + e~ spectral functions. Even though the corrections to the CMD-2 results 
have reduced the discrepancy between (7) and the world average of the direct £>(t~ — > 
v t 7t~7t°) measurements (c/. Section 3.2) from 4.6 to 2.9 standard deviations (adding all 
errors in quadrature), the remaining difference of (—0.94 ± 0.10 r ± 0.26 cc ± 0.11 ra d ± 
0.12su(2)(=t0.32 tota i))% is still problematic. Since the disagreement between e + e~ and r 
spectral functions is more pronounced at energies above 850 MeV, we expect a smaller 
discrepancy in the calculation of a^ ad ' LO because of the steeply falling function K(s). More 
information on the comparison is displayed in Fig. 3 where it is clear that ALEPH, CLEO, 
L3 and OPAL all separately, but with different significance, disagree with the e + e~-based 
CVC result. 

5 Results 

The integration procedure and the specific contributions - near nn threshold, the cu and 
resonances, the narrow quarkonia and the high energy QCD prediction - are treated 
as in our previous analysis [10]. The contributions from the different processes in their 
indicated energy ranges are listed in Table 1. Wherever relevant, the two e + e~- and 
r-based evaluations are given. The discrepancies among them discussed above are now 
expressed in terms of a^ ad ' LO , giving smaller estimates for the e + e~-based data set by 
(-11.9±6.4 cxp ±2.4 rad ±2.6 su(2) (±7.3 total )) lO" 10 for the tttt channel and (-2.8 ±2.6 cxp ± 
0.3 rad ± l-0 S u(2) (±2.9 to tai)) 10~ 10 for the sum of the Air channels. 
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had.LO (10-10) 

Modes Energy [GeV] &+e _ * ^ +e _ 



Low s exp. 7r~^7r — 


[2?7i 7r ± — 


0.500] 


58.04 ± 




[0.500 - 


1.800] 


450.16 ± 


tt^^v try 


[0.500 - 


1.800] 


0.93 ± 


(jj 


[0.300 - 


0.810] 


37.96 ± 


TT^TT~TT® [below </>] 


[0.810 - 


1.000] 


4.20 ± 




[1.000 - 


1.055] 


35.71 ± 


tt^tt~ 7r^ [above (j)] 


[1.055 - 


1.800] 


2.45 ± 


7T + 7T _ 27r 


[1.020 - 


1.800] 


16.76 ± 


27T+27T - 


[0.800 - 


1.800] 


14.21 ± 


27r + 27r _ 7r 


[1.019 - 


1.800] 


2.09 ± 




[1.019 - 


1.800] 


1.29 ± 


37r + 37r _ 


[1.350 - 


1.800] 


0.10 ± 


2tt + 2tt-2tt 


[1.350 - 


1.800] 


1.41 ± 


7r+7r _ 47r ( 3 ) 


[1.350 - 


1.800] 


0.06 ± 


'/ V M / / " 1 


[1.075 - 


1.800] 


0.54 ± 


7r°7)7r° 


[0.975 - 


1.800] 


0.63 ± 


7r° 7 )(7r7r) 


[1.340 - 


1.800] 


0.08 ± 




[1.055 - 


1.800] 


4.63 ± 




[1.097- 


1.800] 


0.94 ± 


tfOtfi^T (3) 


[1.340 - 


1.800] 


1.84 ± 


KKtt (3) 


[1.440 - 


1.800] 


0.60 ± 




[1.441 - 


1.800] 


2.22 ± 


R = J2 excl. modes 


[1.800 - 


2.000] 


8.20 ± 


i? [Data] 


[2.000 - 


3.700] 


26.70 ± 


j/y 


[3.088 - 


3.106] 


5.94 ± 


^(25) 


[3.658 - 


3.714] 


1.50 ± 


i? [Data] 


[3.700 - 


5.000] 


7.22 ± 


-^udsc [QCD] 


[5.000 - 


9.300] 


6.87 ± 


i^dscfc [QCD] 


[9.300 - 


12.00] 


1.21 ± 


Rudscbt [QCD] 


[12.0- 


- oo] 


1.80 ± 



1 70 ± 1 17 56 03 ± 1 60 ± 28 

X • 1 U — 1 — J- • -L 1 ^j KJ , KJKJ — 1 — X • KJKJ — 1 — \J • £jKJ 


+2 ± 2 6 

| ^ • KJ 1 ^ • KJ 


4 89 ± 1 57 464 03 ± 2 95 ± 2 34 


— 13.9 ± 6.4 






1 02 ± 31 




40 ± 05 




84 ± 20 

v7 . Ot: _1_ KJ . jLjKJ 




26 ± 03 

\J • jLj\J — 1 — KJ >KJ*J 




1 31 ± 20 21 45 ± 1 33 ± 60 


-4.7 ± 1.8 


87 ± 23 12 35 ± 96 ± 40 

U«U 1 — 1 — KJ > £JKJ -L — . 'J J — 1 — KJ • iJ KJ — 1 — KJ • i V/ 


+1.9 ± 2.0 


43 ± 04 




22 ± 02 




0.10 ±0.00 




0.30 ± 0.03 




06 ± 00 

V t • V / KJ — 1 — V 7 • KJ KJ 




07 ± 01 

KJ • KJ f — 1 — KJ • KJ -L- 




10 ± 01 




01 ± 00 

V t • V t A — 1 — V / ■ \J KJ 




40 ± 06 

v / • -i v / — i — v / ■ v y ya 




0.10 ±0.01 




24 ± 02 

KJ • ^J^t 1 KJ • KJ Aj 




u.zu ± U.Ul 




1.02 ±0.03 




0.66 ±0.10 




1.70 ±0.03 




0.35 ±0.00 




0.14 ±0.00 




0.28 ±0.00 




0.10 ±0.00 




0.05 ± 0.00 




0.01 ± 0.00 





^ / -i- _ , n x m i 696.3 ± 6.2 CXD 711.0±5.0 exD 1 ._ . _ _ 

^( e + e ^hadrons) 2m w ± - oo + q« +ns Toe -l4.7±7.9 tot 

± 3. b ra d ± U.8 ra d ± 2.8 SU (2) 

1 e + e _ data are used above 1.6 GeV (see Ref. [10]). 

2 Not including u> and <f> resonances (see Ref. [10]). 

3 Using isospin relations (see Ref. [10]). 

Table 1: Summary of the a h ad > LO contributions from e + e~ annihilation and r decays. The 
uncertainties on the vacuum polarization and FSR corrections are given as second errors 
in the individual e + e~ contributions, while those from isospin breaking are similarly given 
for the r contributions. These 'theoretical' uncertainties are correlated among all channels, 
except in the case of isospin breaking which shows little correlation between the 2n and 
47T channels. The errors given for the sums in the last line are from the experiment, the 
missing radiative corrections in e + e~ and, in addition for r , SU(2) breaking. 
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1 

1 

h- 


-•- 
-•- 
-•- 


CLEO 
1 25.42±0. 12±0.42 

OPAL 

1 25.44±0.17±0.29 

L3 

H 25.44±0.16±0.10 

ALEPH preliminary 
1 25.47±0.10±0.09 


eV CVC 

24.52 ±0.32 

, , , 


1 • 

, , 1 




x Average 

25.46±0.10 
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B(x _ ^v x 7t"7t°) (in%) 



Figure 3: The measured branching ratios for t~ — > z/ T 7r~7r compared to the prediction 
from the e + e~ — > 7r + 7r~ spectral function applying the is o spin-breaking correction factors 
discussed in Ref. [10]. The measured branching ratios are from ALEPH [12], CLEO [34] 
and OPAL [35]. The L3 and OPAL results are obtained from their hn branching ratio, 
reduced by the small Kir contribution measured by ALEPH [32] and CLEO [33]. 



The total discrepancy (—14.7 ± 6.9 exp ± 2.7 ra d ± 2.8su(2) (±7.9 to tai)) 1CT 10 amounts to 
1.9 standard deviations. The difference within errors could now be considered to be 
acceptable, however the systematic disagreement between the e + e~ and r tttt spectral 
functions at high energies precludes one from performing a straightforward combination 
of the two evaluations. 



5.1 Results for a M 

The results for the lowest order hadronic contribution are 



a had,LO = (696.3 ± 6.2 cxp ± 3.6 rad ) lO" 10 [e + e~ -based] , 

a had, LO = (711.0 ± 5.0 cxp ± 0.8 rad ± 2.8 SU (2)) 10- 10 [r-based] . 
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Adding to these the QED, higher-order hadronic, light-by-light scattering and weak con- 
tributions as given in Section 2, we obtain the SM predictions 

aJ M = (11 659 180.9 ± 7.2 had ± 3.5 LBL ± 0.4 QED+EW ) 10- 10 [e+e~ -based] , 
af 1 = (11 659 195.6 ± 5.8 had ± 3.5 LBL ± 0.4 QED+EW ) 1(T 10 [r-based] . 

(9) 

These values can be compared to the present measurement (2). Adding experimental and 
theoretical errors in quadrature, the differences between measured and computed values 
are found to be 

a ex P _ a SM = (22.1±7.2 hadiLO ±3.5 othcr ±8.0 exp ) 10- 10 [e + e~ -based] , 

a ex P _ a SM = ( 7 4 ± 5 . 8had LO ± 3. 5othcr ± 8 .o exp ) io- 10 [r-based] , ' ° '' 

where the first error quoted is specific to each approach, the second is due to contributions 
other than hadronic vacuum polarization, and the third is the BNL g-2 experimental error. 
The last two errors are identical in both evaluations. The differences (10) correspond to 
1.9 and 0.7 standard deviations, respectively. A graphical comparison of the results (9) 
with the experimental value is given in Fig. 4. Also shown are our estimates [1,9], obtained 
before the CMD-2 and the new r data were available (see discussion below), and the e + e~- 
based evaluations of Refs. [10,15], obtained with the previously published, uncorrected 
CMD-2 data [11]. 
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Figure 4: Comparison of the results (9) with the BNL measurement [16]. Also given are 
our estimates [1,9] obtained before the CMD-2 data were available. For completeness, we 
show as triangles with dotted error bars the e + e~ -based results [10,15] derived with the 
previously published CMD-2 data [11]. 
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6 Discussion 



Although the new corrected CMD-2 7r + 7r~ results are now consistent with r data for the 
energy region below 850 MeV, the remaining discrepancy for larger energies is unexplained 
at present. Hence, one could question the validity of either e + e~ data with their large 
radiative corrections, r data, or the isospin-breaking corrections applied to r data. We 
shall briefly discuss these points below. 

• The CMD-2 experiment is still the only one claiming systematic accuracies well 
below 1%. It is thus difficult to confront their data with results from other exper- 
iments. Whereas the measurements from OLYA are systematically lower than the 
new CMD-2 results in the peak region, there is a trend towards agreement above, 
as seen in Fig. 2. This behaviour appears to be confirmed by preliminary data from 
the KLOE experiment at Frascati using the radiative return method from the 
resonance [43]. We are looking forward to the final precise results from KLOE and 
from a similar analysis performed by the BABAR Collaboration under very different 
kinematic conditions [44]. 

The relative disagreement between older e + e~ results and CMD-2 can be quantified 
using the CVC prediction: indeed the value of (24.52 ± 0.26 cxp )% obtained for 
#cvc(r _ -> z/ T 7r-7r°), reduces to (23.69 ± 0.68 cxp )% if the CMD-2 data are left 
out, increasing the relative difference with the measured value in r decays from 
(3.8 ± 1.3)% to (7.4 ±2.9)%, a discrepancy hardly compatible with electromagnetic 
isospin breaking. Although the e + e~ data are consistent with respect to the a M 
estimate within their systematic uncertainties, there is some evidence that the older 
data are pulling the value down. 

• The most precise results on the r irir spectral function come from the ALEPH and 
CLEO experiments, operating in completely different physical environments. On the 
one hand, the main uncertainty in CLEO originates from the knowledge of the rela- 
tively low selection efficiency, a consequence of the large non-r hadronic background, 
while the mass spectrum is measured with little distortion and good resolution. On 
the other hand, ALEPH has both large efficiency and small background, the main 
uncertainty coming from the 7T° reconstruction close to the charged pion, necessi- 
tating to unfold the measured spectrum from detector resolution and acceptance 
effects. A comparison of the r spectral functions from ALEPH, CLEO and OPAL is 
given in Fig. 5. Agreement is observed within quoted errors, in particular in the high 
mass region, although CLEO results are a bit closer to e + e _ data there. Overall, 
the r data appear to be consistent. 

• The last point concerns isospin corrections applied to the r spectral functions. The 
basic components entering SU(2) breaking are well identified. The long-distance 
radiative corrections and the quantitative effect of loops have been addressed by 
the analysis of Ref. [14] showing that the effects are small. The overall effect of the 
isospin-breaking corrections (including FSR) applied to the r nn data, expressed in 
relative terms, is (—1.8 ± 0.5)%. Its largest contribution (—2.3%) stems from the 
uncontroversial short-distance electroweak correction [45]. One could question the 
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Figure 5: Relative comparison of the tt + tt~ spectral functions extracted from r data from 
different experiments, expressed as a ratio to the average r spectral function. The lower 
figure emphasizes the p region. 
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validity of the chiral model used. The authors of Ref. [14] argue that the corrections 
are insensitive to the details of their model and essentially depend only on the shape 
of the pion form factor. As the latter is known from experiment to adequate accu- 
racy, it seems difficult to find room for a ~ 10% effect as observed experimentally. 
Nevertheless, considering the situation regarding the first two experimental points, 
it would seem worthwhile to invest more theoretical work into the problem of isospin 
breaking. 

7 Conclusions 

An update of our analysis of the lowest-order hadronic vacuum polarization contribution 
to the muon anomalous magnetic moment has been performed following a reevaluation by 
the CMD-2 Collaboration of their e + e~ annihilation cross sections. Part of the previous 
discrepancy between the e + e~ and r nn spectral functions has now disappeared so that 
the corresponding evaluations of the lowest-order hadronic polarization contribution to 
the muon magnetic anomaly are closer. However, incompatible cross section measure- 
ments remain between 0.85 and 1 GeV so that we do not proceed with an average of the 
two evaluations. The e + e~- and r-based predictions are respectively 1.9 and 0.7 standard 
deviations below the direct measurement from the g-2 Collaboration at BNL. The forth- 
coming results from radiative return with KLOE and BABAR will be decisive to sort out 
the remaining problems in the and An spectral functions. 
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